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The i?-symmetry formalism is applied for the supersymmetric economical SU(3)c® 
SU(3)l <S> U(1)x (3-3-1) model. The generalization of the minimal supersymmetric 
standard model relation among i?-parity, spin and matter parity is derived, and discrete 
symmetries for the proton stability in this model are imposed. We show that in such 
a case it is able to give leptons masses at just the tree level. A simple mechanism for 
the mass generation of the neutrinos is explored. With the new i?-parity, the neutral 
fermions get mass matrix with two distinct sectors: one light which is identified with 
neutrino mass matrix, another heavy one which is identified with neutralinos one. 
The similar situation exists in the charged fermion sector. Some phenomenological 
consequences such as proton stability, neutrinoless double beta decays are discussed. 

PACS. ll.30.Er, 14.60.Pq, 14.60.-z, 12.60. Jv 

1 Introduction 

Although the Standard Model (SM) gives very good results in explaining the observed prop- 
erties of the charged fermions, it is unlikely to be the ultimate theory. It maintains the mass- 
lessness of the neutrinos to all orders in perturbation theory, and even after non-pertubative 
effects are included. The recent groundbreaking discovery of nonzero neutrino masses and 
oscillations [1] has put massive neutrinos as one of evidences on physics beyond the SM. 

The Super-Kamiokande experiments on the atmospheric neutrino oscillations have indi- 
cated to the difference of the squared masses and the mixing angle with fair accuracy [2, 3] 
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Since the data provide only the information about the differences in m^, the neutrino mass 
pattern can be either almost degenerate or hierarchical. Among the hierarchical possibilities, 
there are two types of normal and inverted hierarchies. In the literature, most of the cases 
explore normal hierarchical one in each. In this paper, we will mention on a supersymmetric 
model which naturally gives rise to three pseudo-Dirac neutrinos with an inverted hierarchical 
mass pattern. 

The gauge symmetry of the SM as well as those of many extensional models by themselves 
fix only the gauge bosons. The fermions and Higgs contents have to be chosen somewhat 
arbitrarily. In the SM, these choices are made in such a way that the neutrinos are massless 
as mentioned. However, there are other choices based on the SM symmetry that neutrinos 
become massive. We know these from the popular seesaw [4] and radiative [5] models. 
Particularly, the models based on the SU(3)c <8> SU(3)l ® XJ(l)x gauge unification group 
[6, 7, 8], called 3-3-1 models, give more stricter fermion contents. Indeed, only three fermion 
generations are acquired as a result of the anomaly cancelation and the condition of QCD 
asymptotic freedom. The arbitrariness in this case are only behind which SM singlets put in 
the bottoms of the lepton triplets? In some scenarios, exotic leptons may exist in the singlets. 
Result of this is quite similar the case of the SM neutrinos. As a fact, the mechanisms of 
the Zee's type [5] for neutrino masses arise which been explored in Ref. [9]. 

Forbidding the exotic leptons, there are two main versions of the 3-3-1 models as far as 
minimal lepton sectors is concerned. In one of them [6] the three known left-handed lepton 
components for each generation are associated to three SU(3)l triplets as (vi, /, I c )l, in which 
l c L is related to the right-handed isospin singlet of the charged lepton I in the SM. No extra 
leptons are needed and therefore it calls that a minimal 3-3-1 model. In the variant model [7] 
three SU(3)l lepton triplets are of the form (i/j, I, vf) L , where vf is related to the right-handed 
component of the neutrino field v h thus called a model with the right-handed neutrinos. This 
kind of the 3-3-1 models requires only a more economical Higgs sector for breaking the gauge 
symmetry and generating the fermion masses. It is interesting to note that two Higgs triplets 
of this model have the same U(l)x charges with two neutral components at their top and 
bottom. Allowing these neutral components vacuum expectation values (VEVs) we can 
reduce number of Higgs triplets to be two. Therefore we have a resulting 3-3-1 model with 
two Higgs triplets [10]. As a result, the dynamical symmetry breaking also affects lepton 
number. Hence it follows that the lepton number is also broken spontaneously at a high 
scale of energy. Note that the mentioned model contains very important advantage, namely, 
there is no new parameter, but it contains very simple Higgs sector, therefore the significant 
number of free parameters is reduced. To mark the minimal content of the Higgs sector, this 
version that includes right-handed neutrinos is called the economical 3-3-1 model. 

Among the new gauge bosons in this model, the neutral non-Hermitian bilepton field 
X° may give promising signature in accelerator experiments and may be also the source 
of neutrino oscillations [11]. In the current paper, the neutrinos of the 3-3-1 model with 
right-handed neutrinos is a subject for extended study. 

The 3-3-1 model with right-handed neutrinos gives the tree level neutrino mass spectrum 
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with three Dirac fermions, one massless and two degenerate in mass [12]. This is clearly 
not realistic under the experimental data. However, this pattern may be severely changed 
by quantum effects and gives rise to an inverted hierarchy mass pattern. This is a specific 
feature of the 3-3-1 model with right-handed neutrinos which was considered in Ref.[12] (see 
also [13]), but such effects exist in the very high level of the loop corrections. 

The outline of this work is as follows. In Sec. 2 we define the i?-charge in our model 
in order to get similar results as in the Minimal Supersymmetric Standard Model (MSSM). 
While in Sec. 3 we impose another discrete symmetry that allow neutrino masses but forbid 
the proton decay and the neutron-antineutron oscillation. In Sec. 4 we calculate the fermion 
masses in our model, then we present some phenomenological discussion of this model. Our 
conclusions are found in the last section. In Appendix A, we present the mass matrix 
elements of the neutral fermions. 

2 Discrete R-parity in the supersymmetric economical 
3-3-1 model (SUSYEC0331). 

In the supersymmetric 3-3-1 model with right-handed neutrinos (SUSY331RN) [14], the im- 
parity was already studied and we have shown that if /^-symmetry is broken, the simple 
mechanism for the neutrinos mass can be constructed [15]. This mechanism produces the 
neutrinos mass which is in agreement with the experimental data. 

The supersymmetric extension of the economical 3-3-1 model (SUSYEC0331) was pre- 
sented in [16]. The fermionic content of SUSYEC0331 is the following: the left-handed 
fermions are in the triplets/antitriplets under the SU(3)l group, namely, the usual lep- 
tons are the triplets Li = {yi^vf)^ ~ (1,3,-1/3), i = 1,2,3; while in the quark sec- 
tor, we have two families in the antitriplets Q a i = (d a ,u a , D a ) ~ (3, 3*,0), a — 1, 2, 
and one family in the triplet Q%l = (u 3 ,d3,T) ~ (3,3,1/3). The right-handed com- 
ponents are in the singlets under the SU(3) L group: l\ ~ (1,1,1), u\ ~ (3*, 1,-2/3), 
d\ ~ (3*, 1, 1/3), which are similar to those in the SM. In addition, the exotic quarks trans- 
form as T c ~ (3*, 1, -2/3), D c a ~ (3*, 1, 1/3). 

The scalar content is minimally formed by two Higgs triplets: x — (XiiX~~ iX^Y ~ 
(1,3,-1/3) and p = {pf , p°, pt) T ~ (1)3,2/3). In order to cancel chiral anomalies in 
the SUSYEC0331 model, we have to introduce the followings scalar Higgs triplets x' = 
(X?,X' + ,X2°) T ~ (1,3*, 1/3) and p' = (p'f , p'° , p' 2 ~ ) T ~ (1,3*, -2/3). In this model, the 
VEVs are defined by 

V2( X ) T = (u,0,w), V2( X ') T = (u',0,w'), (2.1) 
V2(p) T = (Q,v,0), V2(p') T = (0,v',0). (2.2) 

The VEVs w and w' are responsible for the first step of the symmetry breaking, while u, vl 
and v, v' are responsible for the second one. Therefore, they have to satisfy the constraints: 

u, u' , v, v' w, w' . (2.3) 
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The complete set of fields and the full lagrangian of SUSYEC0331 are given in Ref. [16]. 
The most general superpotential is given by: 

Wo W-* 

«'- = f + f, (2.4) 



where 



and 



W 2 = PqiUlX + P X XX.' + P P PP \ (2-5) 



W 3 = iabL a Lp'i c bL + \ a eL aL xp + \' ab eL aL L bL p 
+ KiQ 1L x'u c iL + k'Qilx'u'E + ViQiLpd-L 

+ KQlLpd'aL + KaiQaLpUiL + ^'aQaLpuf 

+ ^-aiQaLXdll + R'apQaLXdpL + a^Q (xhQ $lQ^L 

+ iiil3jd c iL d'p L u c jL + t,2i/3d c iL dp L u[ + izijkd\ L d] L u c kL 
+ ^ijd c iL d c jL v!l + ^ a ind'^ L d'p L u c iL + Z,6apd'a L dp L u'l 

+ iaajLaLQaLd C j L + ^ aa pL aL Q aL d!p L . (2.6) 

The coefficients poa,P P and fj, x have mass dimension and can be complex variables [17], 
while all coefficients in W 3 are dimensionless, and \' ab = —X' ba . 

Let us impose the i?-parity as the same as that of the minimal supersymmetric standard 
model. In this case, we have to choose the following /^-charges 

n x = n x' = n p = n p' = O5 
n L = n Qa = uq 3 = 1/2, 

m = n u = n d = n T = n D = -1/2. (2.7) 
The superpotential satisfying the above i?-parity conservation is written as 

W RC = ^XX' + yPP' + I [ X uAp'i C j + mAifui + <Qsv'T c + ^Qsx'ui + k' 2 Q3X'T c 

+ K3aiQ a fjdl + KqapQafjDp + ^iaiQapul + K,' 4a Q a pT C + K 5i Q 3 p ' d\ + K,' 5/3 Q 3 p' Dp 

+ K6aiQaXd C i + KLpQaXDp] ■ ( 2 -8) 



With this superpotential, we have shown that [16] the boson, Higgs sectors and the fermion 
one gain masses. 

Thus, the i?-parity in this model, as in the SUSY331RN, can be re-expressed via the spin 
S, new charges C and B in terms of [15] 

J R-parity= (-l) 25 (-l) 3(B+/:) , (2.9) 
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where the charges B and C for the multiplets are defined as follows [12] 
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(2.11) 
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From the superpotential given in Eq.(2.8), it is easy to see that the charged leptons gain 
mass through the term 

Xuj 



Up I] + he. 



(2.13) 



Their mass matrix, see Eq.(2.2), is given by 



V A 



(2.14) 



Note that there is only VEV of p' given the charged leptons masses. 

Unfortunately, as in the MSSM case, due to conservation of the i?-parity defined in 
Eq.(2.7), there are no term which gives neutrinos masses. However, looking at the superpo- 
tential of this model given in Eq.(2.6), we see that there is a term LiLjp, which generates 
the following 

(2.15) 



(2.16) 



\ 3ij [LiLjp + LiLjp + LiLjp 
The first term in (2.15) generates the following neutrino mass matrix 

V A 



As shown in Ref. [16], the mass pattern of this sector is 0, 0, m u , m u , m u , m u . Note that in 
this case, we have two massless neutrinos. Unfortunately, as in the non-symmetric version, 
the quantum corrections at one loop level cannot generate the realistic mass spectrum to 
the neutrinos. To get the realistic neutrino masses, one have to introduce new physics scale 
or inflaton with mass around the GUT scale [18]. 

In this article, we will explore a new mass mechanism to generate neutrino masses at tree 
level for all neutrinos and study the flavor violating processes, such as neutrinoless double 
beta decay, which do not exist in our previous work. 
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3 The discrete symmetry for proton stability and neu- 
trino masses in SUSYEC0331 



In this section to get neutrino mass and impose flavor violating processes, we chose the new 
R charge as follows 



n L = n x 



n x > = n u = tit 



2' 

1 

~2' 

n Q A = u p' = !» n d = n D = -2, 
n p = -1, n Qa = -, m = --. (3.1) 

This i?-charge is different from those presented in Ref. [16]. We will show that in this case, 
there exist some new phenomena, which are previously not allowed. 

The terms under this symmetry are obtained by 

W = W RC + ^oiLiX + - (A 2i eLiX/3 + X^eLiLjp + X'^jQaUdj + ^aipQaUDp) , 

(3.2) 

where Wrc is defined in Eq.(2.8). The superpotential given in (3.2) will not only allow 
some interesting flavor violating processes but also will simultaneously give the nucleons a 
stability. As we will show in next section, this superpotential will also generate masses to 
all neutrinos in the model. 



4 Fermion masses 

The superpotential (3.2) provides us the mixing between the leptons and higgsinos as 

-tfLtf - ^f{Lap + pU X ) + H.c. (4.1) 

With the above terms, we get the mass matrices for the neutral and charged fermions. 
Diagonalizing these matrices we obtain the physical masses for the fermions. Firstly, let us 
study the neutral fermions masses. 

4.1 Masses of the neutral fermions 

In the basis \&° of the form 

(u, v 2 v z v\ vl vlx\ Xi Xl X' 2 ° ~P° P'° BmmX X*) , 



the mass Lagrangian can be written as follows 

--^°) T Y ^ + H.c. (4.2) 

Here Y° is symmetric matrix with the nonzero elements given in Appendix.A, where the 
mass eigenstates are given by 

X? = N^, .7 = 1,..., 17. (4.3) 

The mass matrix of the neutral fermions consists of three parts: (a) The first part M v 
is the 6x6 mass matrix of the neutrinos which belongs to the SUSYEC0331; (b) The 
second part is the 11 x 11 mass matrix of the neutralinos, which exists only in the 
presented supersymmetric version, has been analyzed in [19]; (c) The last part M v n arises 
due to mixing among the neutrinos and the neutral higgsinos. Thus, the mass matrix for 
the neutral fermions is signified as follows 

Y = 1 KU*n <"«>..... j ' <4 ' 4) 

where matrices M v and M vN are presented in Eq.(A.l) and Eq. (A. 3). 

Let us keep the mass constraints from astrophysics and cosmology [20] as well as being 
consistent with all the earlier analysis [21], the parameters in the mass matrix can be 
chosen as a typical example: 

Hp = 600 GeV, fi x = 700 GeV, 
M 3 = M 8 = 300GeV, M 45 = 400 GeV, 
A*oi = At02 = H03 = 1 GeV, 

A2I = A22 = A03 = 1, 

A312 = 4-KT 11 , A 3 i 3 = 5-10- n , 
A321 = 6-KT 11 , A 323 = 7- 10- 11 , 

A331 = 8-10~ n , A332 = 9 • 10- 11 , (4.5) 

Here in this model, the Higgs bosons' VEVs are fixed as follows 

v Xl = 15GeV, v x[ = 10 GeV, 
v p = 244.9 GeV, v p , = 13 GeV, 

v X2 = Vx , 2 = 1000 GeV, (4.6) 

and the value of g is given in Ref. [20]. 

Using the values given in Eqs. (4.5,4.6), the eigenvalues of fermion mass matrix are ob- 
tained as 

77i*o = -1282,27 GeV, m^o = 1236,53 GeV, m^o = 1041 GeV, m*o = -705,43 GeV, 

X17 ' Xl6 ' X15 ' X14 
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77V) = 631,46 GeV, m*o = 620, 21 GeV, rn^o = 487, 49 GeV, m^o = 385, 08 GeV, 

m^o = 304, 82 GeV, m^o = 186, 45 GeV, m^o = 57, 13 GeV, m^o = -0, 103 GeV, 

m .o = -0, 043 GeV, m^o = 2, 0415 • 10~ n GeV, m^o = 2, 0413 • 10~ n GeV, 

777^0 = -2, 0412 • 10~ n GeV, m^o = -2, 0410 • 10~ n GeV. (4.7) 

In the Eq.(4.7), there are some negative eigenvalues. In order to obtain the positive mass, 
the eigenstates need to be redefined by the chiral rotations. 

Eq.(4.7) shows that we have two very distinct sector, one contains the light neutral 
fermions that we will associate with the usual neutrinos in the SM and the other one contains 
the heavy neutralinos. The lightest neutralino mass equals to 57 GeV and it is consistent 
with limits on inelastic dark matter from ZEPLIN-III [22]. 

Using the values given in Eqs. (4.5,4.6), the eigenvalues of the mass matrix (M„) 6x6 are 
obtained as 

777*0 = 777^0 = GeV, 

Xi X2 

m -o = 777*0 = 777^0 = 777*0 = 2, 125 • 10 -20 GeV. (4.8) 

These values are smaller than that of the new mechanism given in (4.7). On the other 
hand, the eigenvalues of the matrix (Mrv) llxll are obtained by putting the numerical given 
in Eqs. (4.5,4.6) as follows 

777*0 = 1207 GeV, m*o = 1143 GeV, m*o = 1040 GeV, m*o = 704, 61 GeV, 

X7 ' Xg 7 Xg 7 Xio 7 7 

m*o = 585,49 GeV, m*o = 499, 85 GeV, m*o = 429, 68 GeV, m*o = 374, 72 GeV, 

All ' Xl2 Xl3 ' X14 

m - ?s = 304,81 GeV, m*^ = 175,89 GeV, m*^ = 56,88 GeV. (4.9) 

These results can be understood as follows: Because of the interference matrix (M vN )q x u 
between the neutrino mass matrix and the neutralino mass matrix, all neutrinos gain mass 
at the tree level. This change of the neutrino mass spectrum is suitable to experiment data. 
Thus the neutrino mass spectrum in the model under consideration depends on the choice 
of .R-parity. Now we deal with the charged fermions. 

4.2 Masses of the charged fermions 

To write mass matrix of the charged fermions, we will choose the following bases 

r = ( h k h w- y- p'f fc r ) T , 

^ = (ii i% it w + y+ pf pt f + ) T , (4.io) 

and define 

#± = ( V + ^~ Y ■ (4.11) 
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With these definitions, the mass term is written in the form 



1 



— (^ ± ) J F ± ^ ± + H.c, 
2 
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(4.12) 



(4.13) 



The chargino mass matrix Y ± is diagonalized by using two unitary matrices, D and E, 
defined by 



X 



+ 



= DijV+ Xi =E ij Vj, i,j = 1, 



(4.14) 



The characteristic equation for the matrix Y ± is 



det(F ± - XI) = det 

Since X T ■ X is a symmetric matrix, A 2 must be real and positive because Y ± is also sym- 
metric. In order to obtain eigenvalues, one only have to calculate X T ■ X. The diagonal 
mass matrix can be written as 

M SCM = E*XD- 1 . (4.16) 
To determine E and D, it is useful the following observation 

M 2 SCM = DX T ■ XD- 1 = E*X ■ X T (E*)-\ (4.17) 

It means that D diagonalizes X T ■ X, while E* diagonalizes X ■ X T . In this case we can 
define the following Dirac spinors: 

= ( xt 11 f, * c (xr) = ( Xi It ) T ■ (4.18) 
where xt ls the particle and Xi is the anti-particle [23]. 



-A X T 
X -A 



det(A 2 -X T -X). 



(4.15) 
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Using the values given in Eqs. (4.5,4.6), the eigenvalues of the charged fermion matrix 
given at Eq.(4.14) are obtained as 

m e = 5 x KT 4 GeV, = 0, 104 GeV, m r = 1, 179 GeV, 
77i ± = 69, 09 GeV, m-± = 522, 92 GeV, 

m ± = 600,22 GeV, m-± = 732,64 GeV, m~± = 1168 GeV. (4.19) 

A3 A4 X5 v ' 

It is easily to see that the mass matrix is divided in two sectors: one heavy which is 
identified as the charginos and has been studied in our previous work [16]. Another one light 
which is identified as the usual leptons. 

On the other hand, if we take the mass matrix from our work and using the same values 
to the parameters we get masses of the charged leptons: 

m e = 5 x 10~ 4 GeV, m M = 0, 119 GeV, m T = 1, 249 GeV. (4.20) 

Now using the mass matrix to the charginos given in [19] we get 

m - i = 69, 05 GeV, m^ = 522, 78 GeV, m^ 3 = 600, 01 GeV, 

m i4 = 600,42 GeV, m* 8 = 1168 GeV.GeV. (4.21) 

In this case, the new elements put down masses of the muon and tauon and at the same 
time remove the mass degeneracy between X3 an d Xi- Anyway we can see that the mass 
matrix in the charged fermion sector is basically divided in two sectors: one giving masses 
of the usual known leptons and another one giving masses of the new charginos. 



5 Nucleon decay in the SUSYEC0331 

Let us remind that in the MSSM [24, 25, 26, 27, 17], the i?-parity violating terms in super- 
potential are given by W 2 rv + W 3 rv + W 2RV + W 3RV) where 

W 2 RV = HoatL a H 2 , 

W 3RV = X abc eL a L b l c c + yl iai eQiL a d) + A-^d-it^. (5.1) 

Here we have suppressed SU(2) indices, e is the antisymmetric SU(2) tensor. Above, and 
below in the following, the subindices a, b, c run over the lepton generations e, /z, r but a 
superscript c indicates charge conjugation; i,j, k = 1,2,3 denote quark generations. 

Note that the A'-coupling in the MSSM is similar to the ^-coupling in the most general 
form of the superpotential in the SUSYEC0331, which is given in Eq.(2.6), and the coupling 
A" is similar to the coupling £ 3 . 

From Eq. (5.1), we can obtain the B- violating Yukawa couplings as follows 

\'; jk d lR u c jL d c k + H.c. (5.2) 

10 



The interactions among lepton, quark and squark are given by 

Kaj{d c iR v aL -u c iR l aL )d j + H.c. . (5.3) 

Taking into account Eqs. (5.2,5.3), we can draw two Feynmann diagrams describing the 
proton decay, which are shown in the Fig.l and Fig. 2. The Fig.l describes the proton decay 
into charged leptons. At the lowest approximation, there is no mixing in the quark, neutrino 
and squark sectors. It means that u± = u, ui = c and u% = t and so on. The proton could 
decay into p — > n°e + , p — > D°fi + and p — > (■u 1 t c )r + , but the last two modes are forbidden 
kinematically. 




Figure 1: Proton decay into charged leptons in the MSSM and in SUSYEC0331 (with 
A" ^£3 and A'->&). ' 

The analysis presented above shows that the proton can decay only in p — > n°e + . On 
dimensional grounds, we estimate 

r(p ttV) « a{Klk)a } Klk) Ml oton) (5.4) 

where a(X) = A 2 /(4tt). Giving r(p en) > 1.6 x 10 33 years [20] and taking m~ dk ~ C(lTeV), 
we obtain 

A' nfc A? lfc < 5.29 x 10- 26 . (5.5) 

It is consistent with the limits presented in [28]. For a more detailed calculation see [29, 30]. 
Other decay modes, where the proton decay into antineutrino, have been considered in [31]. 
The mentioned decay modes are p — > n + i7 e , p — y K + i7^, p — > B + i7 T 1 . In these cases, we get 
the same numerical results as presented in Eq.(5.5). 



1 Again, we are bare mixing in the quarks, neutrinos and squarks sectors 
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Figure 2: Proton decay into antineutrinos in the MSSM and in SUSYEC0331 (with A" — > £3 
and A' ->■ £7). 

The bound presented in Eq.(5.5) is so strict. So the natural explanation only is that at 
least one of the couplings has to be zero. In order to avoid that the simplest way is to impose 
the i?-symmetry. This leads to avoid the proton decay. 

In our superpotential given at Eq.(3.2), we allow only interactions that violate L-number. 
Therefore the proton is stable at tree-level. However it is not only forbidden the danger- 
ous processes of proton decay but also forbidden the neutron-antineutron oscillation. This 
oscillation was studied in detail in [32, 33, 34, 35]. 

In [14, 15], in the framework of the supersymmetric 3-3-1 model with right handed 
neutrinos, the R-parity violating interaction was applied for instability of the nucleon. The 
result is consistent with that of the present article (noting that in Ref [14] , the authors have 
taken a lower limit of the proton lifetime equal to 10 32 years ). 



6 Neutrinoless double beta decay in SUSYEC0331 



Neutrinoless double beta (Ov/3/3) decay is a sensitive probe of physics beyond the stan- 
dard model, see Fig. (3), since it violates the conservation of lepton number [36, 37, 38, 39]. 
Such experimental results if available will cast stringent constraints on new physics. The 
nucleon level process of 0vf3f3 decay, n + n — >■ p + p + e~ + e~, can be obtained via the 
lepton number violating subprocess, d + d — >u + u + e + e. For experiments, the parent 
nucleus Az can decay to a daughter nucleus Az+2 via the two steps of virtual transitions, 
Az — > Az+i + j3~ — > Az+2 + P~ + j3~ . It is to be noted that only the ordinary double 
beta reactions have been experimentally observed, while the active searches for the Qv/3/3 
reactions are actually get with null results. The experimental measurements of 0vf3f3 de- 
cay are currently accounted for even-even heavy nuclei such as 48 Ca — > 48 Ti, 76 Ge — > 
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Figure 3: Neutrino less double beta decay in SM with massive neutrinos. 

76 Se, 82 Se ->■ 82 Kr, 100 Mo ->■ 100 Ru, 128 Te -> 128 Xe. The most stringent bounds have 
been given by the Moscow-Heidelberg collaboration [37, 38] on the Ov/3/3 decay half-life of 
76 Ge as T ud > [1.1 x 10 25 , 1.5 x 10 25 ] yr, respectively. We can find in Ref. [40] for a general 
review, and Ref. [39] for a summary of the experimental results as well as the future projects. 

Particularly, for the conventional 0z//3/3-decay with massive Majorana neutrino exchange 
it implies an upper bound on the neutrino mass below 1 eV. The supersymmetric mechanism 
of 0vj3j3 decay was first suggested by Mohapatra [41] and further studied in Refs. [42, 43]. In 
Ref. [44] , the R parity violating Yukawa coupling of the first generation is strongly bounded 
by A' m < 3.9 ■ 10" 4 due to the gluino exchange (V/3/3-decay. Babu and Mohapatra [45] have 
latter implemented another contribution comparable with that via the gluino exchange. This 
set stringent bounds on the products of R parity violating Yukawa couplings A' 1:li A' ia of ith 
generation index, see Fig. (4). The constraints on A' m coming from the Ou/3/3 half-life limit 
have been calculated [46, 47] 

with the assumption m d ~ R ~ m Si and m,g, rrig being squark and gluino masses, respectively. 
We find further [46] 

An 3 A' 13 i < 1-1 -10- 7 , (6.2) 
A'mA'm < 3.2 • 1(T 6 . (6.3) 

The analyses presented in the MSSM to the Oz//3/3-decay are still hold in the SUSYEC0331 
model, because the A' is allowed in our superpotential as shown at Eq.(3.2). 
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Figure 4: Neutrinoless double beta (0vj3j3) decay in the MSSM with massive neutrinos and 
in SUSYEC0331 (with A' ->■ £ 7 ). 

7 Conclusion 

In this paper we have presented new i?-symmetry for the supersymmetric economical SU(3)c® 
SU(3)l ® U(l)x mode and studied neutrino mass by implication for the obtained -R-parity. 
The neutrino mass spectrum is affected by the chosen R-parity. By imposing the i?-parity, 
namely: R-parity = (— l) 2S (— l)^ B+c \ the neutrino mass spectrum at the tree level contains 
two massless. 

We remind that in the non-supersymmetric economical 3-3-1 model, to give neutrinos a 
correct mass pattern, we have to introduce new mass scale around the GUT scale. The same 
situation happens in the supersymmetric version and the above puzzle can be solved with 
the help of the inflaton having mass in the range of the GUT scale. In this paper, we have 
showed that by the chosen R parity and the set of parameters, the pseudo-Dirac neutrino 
mass is available. 

We have found that the other .R-parity given in (3.1) leads to the interference between 
the neutrino and neutralino mass matrices. Because of this interference mass matrix, all 
neutrino gain mass only just at the tree level and the interference mass matrix does not 
much affect on the neutralino mass spectrum. In the charged lepton sector, with the new 
R parity, there is also an interference mass matrix between the usual leptons and charginos. 
By taking the numerical, we show that the charged sector is basically divided in two distinct 
sectors: one giving the usual known leptons and the other ones given the new charginos. 
If we ignore the interference charged lepton mass matrix, the chargino mass spectrum is 
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degenerated. This degenerated mass spectrum is removed by imposing the new R-parity. 

The new i?-parity not only provides a simple mechanism for the mass generation of the 
neutrinos but also gives some lepton flavor violating interactions at the tree level. This will 
play some important phenomenology in our model such as the proton's stability, forbiddance 
of the neutron-antineutron oscillation and neutrinoless double beta decay. 
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A Elements of neutrino mass matrix in SUSYEC0331 

The elements of Y° presented in Eq. (4.2) is given by Eq.(4.4) where 
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